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The problem of separation of different signals in the Cosmic Microwave Background 
(CMB) radiation using the difference in their statistics is analyzed. Considering samples 
of sequences which model the CMB as a superposition of signals, we show how the 
Kolmogorov stochasticity parameter acts as a relevant descriptor, either qualitatively or 
quantitatively, to distinguish the statistical properties of the cosmological and secondary 
signals. 

1. Introduction 

In the study of the Cosmic Microwave Background radiation, one of the basic in- 
formation sources on the early Universe, one deals with a principal problem of 
separation of the cosmological signal from perturbations of various nature 21. The 
CMB reflects not only effects occurring at the early or large scale Universe, the 
Sachs- Wolf and Sunyaev-Zeldovich effects, but also more local perturbations of the 
Galactic, interplanetary, and up to instrumental origin. 

Various descriptors are applied to separate the CMB cosmological signal from a 
secondary signal of particular nature. In the present paper we will apply the Kol- 

|9|oi 

mogorov stochasticity parameter - to reveal the difference in the statistical prop- 
erties of contributions irrelevant to the main CMB signal. As was already shown, 
the Kolmogorov's parameter is an informative descriptor when applied to the CMB 
temperature sequence s, n amely, it enables one to locate regions of different random- 
ness in the sky maps^^^. For example, the Cold Spot, the non-Gaussian anomaly 
in the Southern sky, not only shows enhanced value of Kolmogorov's parameter but 
also possess its variation along the rad ius peculiar to the voids in the large scale 
matter distribution in the U niver se the void nature of the Cold Spot is among 
the discussed interpretations l^^^. 

We compute the Kolmogorov stochasticity parameter and then Kolmogorov's 
distribution for real CMB maps, namely, for the 3 bands of the Wilkinson Mi- 
crowave Anisotropy Probe's (WMAP's) temperature data in which the nature and 
the level of the Galactic radiations is different. We also estimate the simulated maps 
which all together indicate that Kolmogorov's parameter indeed is different for maps 
with different contamination. This motivates the interest to model the behavior of 
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Fig. 1. 'I'(A) maps calculated for WMAP W,V,Q-band temperature data (a-c), respectively, and 
for the simulated, no-noise (d) CMB maps. 



Kolmogorov's parameter for the CMB (continuing L!^, via specially constructed se- 
quences composed of subsystems of different statistic, i.e. with expected statistical 
properties similar to the perturbations present in the primary signal in CMB maps. 
Such 'numerical laboratory' results show that in certain cases one can clearly distin- 
guish the signals due to the qualitative differences in the behavior of Kolmogorov's 
parameter, in other cases, the quantitative properties are the indicative ones. 



2. Kolmogorov statistic 

First let us recall the definition of the Kolmogorov stochasticity parameter 
Consider a finite random sequence of real numbers xi,X2, ■ ■ ■ ,Xn, so that are 
sorted in an increasing manner xi < X2 < ■ ■ ■ < Xn- 
The empirical distribution function is defined as 



F„ (X) = {relative normalized number of the elements Xi which are less than X), 
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I 



September 24, 2009 8:46 WSPC/INSTRUCTION FILE kmodeM 



so that 



FniX) = { 



, X <xi ; 

k/n , Xk < X < Xk+i, fc = 1, 2, . . . ,n — 1 

1 , Xn<X. 



The theoretical distribution function F{X) is defined as 

F{X) = {probability of the event X < X). 
Kolmogorov's stochasticity parameter A„ for a sequence of n values of x is 

X„ = ^-sup\F{X)-Fn{X)\, (1) 
so that Xn itself is a random variable. 

According to the Kolmogorov's theorem A„ is a random number having empirical 
distribution function 

$„(A) = {relativenormalizednumber of the elements Xn less than A), 

uniformly converging to a function $(A) at n — > oo for any continuous distribution 
F{X) 

$„(A) ^ $(A), (2) 

where 



+00 

$(A)= J2 (-1)^6-2'='^', -oo<fc<+oo. (3) 

fe=— oo 



$(A) is Kolmogorov's distribution, it varies from $(0) = to $(00) — 1 mono- 
tonically. 

Tabic 1 shows the Kolmogorov's distribution $(A) and its derivative ^(A) = 
$'(A) at the values of interest of A. 
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A 


$(A) 


0(A) 


0.2 


5.05- 10-13 


1.53-10-1'' 


0.4 


0.0028 


0.1012 


0.6 


0.1357 


1.3241 


0.8 


0.4558 


1.627 


1 


0.73 


1.071 


1.2 


0.8877 


0.5385 


1.4 


0.9603 


0.2222 


1.6 


0.988 


0.0764 


1.8 


0.996932 


0.022 


2 


0.999329 


0.0053 


2.2 


0.999875 


0.0011 


2.4 


0.99998 


0.00019 



From Table 1 we see that the probable values of A are between 0.4 and 1.8. 

According to the definition, for large enough n and random sequence Xn the 
Kolmogorov stochasticity parameter A„ will have a distribution close to $(A). If 
the sequence is not random, the distribution will be different. So, <I'(A) denotes the 
degree of randomness of a sequence. 

For each studied sequence one can get only one value of A„. But that is not 
enough to conclude on the degree of randomness of that sequence. Therefore, con- 
sidering sequences of large enough length n, we split them into subsequences, and 
then calculate A„ for each subsequence. 

Therefore by our strategy we will consider the splits of the sequence Xn into m 
subsequences. In each case we deal with new sequences A™ of length m, which should 
be a random sequence of a distribution close to ^(A) according to Kolmogorov's 
theorem. 

3. CMB $ maps 

We computed the Kolmogorov stochasticity parameter A„ for the CMB temper- 
ature maps, then obtained the Kolmogorov's distribution <I> over the sky. CMB 
temperature maps are given in HEALPIX coordinate system that can have differ- 
ent resolutions described by parameter NSide. We have used WMAP's CMB maps 
of NSide — 512, in order to have enough temperature sequences to obtain lower 
resolution NSide = 32 $ maps, i.e. we reconstruct a sequence from pixel s in one 
cell of lower resolution map and calculate A then for more details see We 
used the WMAP's 5-year data-in the W,V and Q bands (Fig.l a-c). It is known 
that the contamination of the cosmological CMB signal is different in these maps 
due to the different cont ribut ion mainly of the Galactic synchrotron and the dust 
radiation (for details see and references therein). 

As it seen from Figure la-c the $-maps are different for different bands, which 
shows that Kolmogorov's distribution is indeed sensitive to the degree of the contam- 
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ination by specific signals. Then, in Fig. Id we show the <I>-map for the temperature 
map simulated by a standard scheme for the power spectrum parameters, i.e. the 
purely cosmological signal, without any noise. The difference of the simulated $- 
map from those of real maps is also visible. This indicates the principal possibility 
of the simulation of a given signal or several signals and then the comparison of 
their <i>-maps with the observed ones, in order to reveal the contribution of a given 
component. However, to do that one has to model signals with various compositions 
of sub-signals, to probe their sensitivity to the Kolmogorov's statistic. 

4. The sequences and their properties 

We have analyzed sequences of uniform distribution: 7000 different sequences of 
length 10000 each were generated, which then were divided into 7 groups, each of 
1000 sequences. 

The list of the considered sequences is as follows: 

I: random; 

II: non-random: 

Xn = 107 n mod 513; (4) 
III: sequences including two copies of the same random sequence, 

Xi,X2, ■ ■ ■ ,Xn/2:^li^2: ■ ■ ■ :Xn/2'j (5) 

IV: four copies of a same random sequence; 
V: random with perturbations, 

= esin^i, (6) 

with parameter e = 0.03; 

VI: the same as V but with e — 0.01; 

VII: non-random, where the two neighbors are equal: 

Xi,Xi,X2,X2,.-.,X„,Xn- (7) 

Figure [2] shows the results of the computations. In all figures the horizontal axis 
is A, the vertical axis denotes $(A), the dashed line is the empirical distribution 
function, the solid line is the $(A). The figures shows the level of coincidence of 
^(A) with the empirical distribution function for each group. 

We see that only for random sequences ^(A) practically coincide with the em- 
pirical distribution function, as it is expected by the Kolmogorov's theorem. In all 
other cases the $(A) deviates from the empirical distribution function. For example, 
for non-random sequences of type Xn — 107 n mod 513 ^(A) is quite different from 
the empirical distribution function, as the values of A vary between 0.195 and 0.230. 
For sequences including two and four copies of the same random sequence, the val- 
ues of A are, respectively, -\/2 and 2 times larger than for those for the random case 
due to the factor i/n in the expression for A. In the perturbed case with e = 0.03 
we have much bigger values of A, so <i>(A) = 1. 
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Table 2 contains ^{X) and <i>„(A) for the I group, the random sequences. 



A 


a>(A) 


$„(A) 


0,38996 


0,00193 


lE-3 


0,5707 


0,09945 


0,101 


0,63823 


0,19 


0,201 


0,69643 


0,28284 


0,301 


0,76231 


0,39352 


0,401 


0,82529 


0,4964 


0,501 


0,88579 


0,58735 


0,601 


0,97301 


0,69994 


0,701 


1,07263 


0,7999 


0,801 


1,23134 


0,90361 


0,901 



Then we split the sequences into 50 subsequences, obtain 50 values of A for 
each subsequence, and calculate Xmean- We are interested in the difference between 
the values of Xmean and theoretical mean XMean = / X(j){X)dX « 0.875029 where 
(^(A) = ^''(A). For random sequences Xmean must be equal to theoretical mean value 
XMean ~ 0.875029. Therefore this difference is a direct measure that describes the 
degree of the randomness o f a se quence, bigger is that difference, more regular is 
the sequence, for details see^^IEl 

In Figures [3] and [4] horizontal axis i is the index of the sequence, so each 1000 
points corresponds to one group noted above. Figure [3] shows that, the difference 
between groups I and III or IV are not obvious, i.e. the split subsequences do not 
"feel" the duplicates. Also, for V and VI groups the perturbations are small, so 
again there is no visible difference between them and the group I. The random case 
differs from II and VII groups, because of the different scales. 

Then we found out the behavior of for each group of sequences, is calcu- 
lated for F{X) and F„(X) and shows how far is the distribution from the uniform 
one. Here is defined as 

^ -Z. • (8) 

i=l 

To generalize the considered above results, we studied sequences of the type 

Zn ^ axn + {I - a)yn, (9) 

where a between 1 and specifies the limiting the random Xn, and regular of 
the form y„ = 107z r7io(i513, cases, respectively. 101 sequences were studied, each 
containing 100 subsequences of length 10000. Steps of 0.01 were used for a, so, to 
each a correspond 100 sequences and hence we have 10100 sequences. Plots of <&(A) 
vs A are shown for different values of a in Figure O 



I 



I 
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From the Figure [Bfa) we can see that when a = (the regular case), the normal- 
ized Xmean precisely corresponds to the normalized value for the second group of 
Figure [3] For a = 1 (the random case), the normalized Xmean precisely corresponds 
to the normalized Xmean for the second group shown in Figure [31 

In Figure [6] and normalized Xmean for a = and a = 1 cases behave the 
same way as in 2nd and 1st groups respectively in figures 2] and [31 

We were also interested in the behavior of Xmean depending on the parameters 
of sequences i.e. the number of sequences, their length, number of the splits, type 
of yn, as well as on the steps of a. 

We also changed the type of their length and the number of splits, then 
fitted the corresponding numerical results with an exponential decaying function. 
Figure [71 shows the results of computations. The left plots are the normalized Xmean 
dependences on the index of the sequence, the right plots are the corresponding 
fits. In each plot the legend above indicates the type of the generated sequences, 
the length and number of splits: for example, the Figure [Tja) is a result of the 
generated sequences of the type 107i mod 513, 10000 is the number of elements in 
each sequence and 50 is the number of splits of the sequences. 

One can see that the maximum of 

\Xmean " 0.875029| (10) 

depends on the parameters of y„, the length of sequences and the number of splits, 
however, in the cases where the length and the number of splits ratios are equal, 
maximal values of \Xmean — 0.875029| are also equal to a precision about lO"-^. 

When larger length sequences were applied, the error bars of the belt in Figure[71 
decrease and due to the large numbers of points the value of increases, as shown 
in the legends. 

5. Conclusions 

We analyzed several types of sequences using the Kolmogorov stochasticity pa- 
rameter to model the contributions of signals of different statistic in the Cosmic 
Microwave Background radiation. The previously undertaken estimations of the 
Kolmogorov's parameter for the real CMB temperature maps obtained by WMAP 
had revealed the potential of the meth od in the separation of the Galactic disk and of 
other structures in the previous works 1^^. We have also computed the Kolmogorov 
CMB maps for various frequency bands of WMAP, which show the variation of 
the statistic at various degree of contamination of the signal by the Galaxy. Those 
peculiarities in the properties of the CMB temperature anisotropics have defined 
the strategy of our present study. Namely we did as follows: once we saw that the 
Galactic contribution can be separated in CMB maps with Kolmogorov analysis 
due to its different degree of randomness, we modeled systems comprised of signals 
that differ in their randomness. Such strategy in principle can help to reveal the 
role of a given sub-signal. 
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Although our numerical experiments deal with specific sequences, the performed 
analysis shows that the CMB can be in principle be separated into sub-signals using 
the Kolmogorov statistic and the stochasticity parameter, either qualitatively or 
quantitatively. 

Some of the systen is con sidered above can serve as numerical illustrations to the 
examples discussed in 13111. 

We are grateful to V.G. Gurzadyan, A. A. Kocharyan, T. Ghahramanyan and G. 
Yegorian for discussions and help and to the referee for many valuable comments. 

References 

1. G. Hinshaw, et al, ApJ Suppl. 180, 225 (2009). 

2. A.N. Kolmogorov, G.Ist.Ital.Attuar, 4, 83 (1933). 

3. V. Arnold, Nonlinearity, 21, T109 (2008). 

4. V.G. Gurzadyan, A. A. Kocharyan, A&A, 492, L33 (2008). 

5. V.G. Gurzadyan, A. A. Kocharyan, A&A, 493, L61 (2009). 

6. V.G. Gurzadyan, et al, A&A, 497, 343 (2009). 

7. V.G. Gurzadyan, et al, A&A, 498, LI (2009). 

8. L Masina. A. Notari, larXi? :0905.1073 

9. G. Rossmanith, et al, larXiv :0905.2854 

10. T.Ghahramanyan, et al, Mod.Phys.Lett. 24, 1187 (2009). 

11. B.T.Draine, A.Lazarian, ApJ, 494, L19 (1998). 

12. V.L Arnold, ICTP/2008/001, Trieste. 



September 24, 2009 8:46 WSPC/INSTRUCTION FILE kmodeM 




I 



I 



September 24, 2009 8:46 WSPC/INSTRUCTION FILE kmodeM 



10 



0,8-1 
0,7- 
0,6- 




' T ' I ' I ' I ' I ' I ' I ' I ' I ' I 

-1000 1000 2000 3000 4000 5000 6000 7000 8000 



Fig. 3. Normalized Xmean for all sequence groups, i denoting the number of the 
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Fig. 5. Kolmogorov's distribution $(A) for Zn = axn + (1 — a)yn for (1) a = 0.2, (2) a = 0.4, 
(3) a = 0.6, (4) a = 0.8. As in Fig. 3, the solid line denotes the Kolmogorov distribution (3), the 
dashed line denotes the computed empirical distribution. 
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